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We develop a semiclassical approximation for the spectral Wigner and Husimi functions in the 
neighbourhood of a classically unstable periodic orbit of chaotic two dimensional maps. The predic- 
tion of hyperbolic fringes for the Wigner function, asymptotic to the stable and unstable manifolds, 
is verified computationally for a (linear) cat map, after the theory is adapted to a discrete phase 
space appropriate to a quantized torus. The characteristic fringe patterns can be distinguished even 
for quasi-energies where the fixed point is not Bohr-quantized. The corresponding Husimi function 
dampens these fringes with a Gaussian envelope centered on the periodic point. Even though the 
hyperbolic structure is then barely perceptible, more periodic points stand out due to the weakened 
CNJ . interference. 
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I. INTRODUCTION 



The mixture of a large number of quantum states within a narrow energy window can be in principle well described 
by semiclassical mechanics |lj|. Yet it is not trivial to predict which classical structures will stand out within the 
interference pattern obtained by the various trajectories that contribute to a given quantum representation of the 
^ \ mixed state. So far, most attention has been directed to the immediate neighbourhood of periodic orbits. For an 
energy such that the action of the orbit is nearly Bohr-quantized, there is an intensity enhancement, known as 
a scar. 

Recently, it has been shown that a periodic orbit in a 4-D phase space actually generates a 2-D ring pattern in the 
Wigner function which may be more clearly distinguishable overall than the scar in its immediate neighbourhood M . 
These rings will be absent in a 2-D map, where the periodic orbit is reduced to a set of isolated points. However, such 
a mapping is the simplest case in which to study the imprints of the stable and unstable manifolds on the Wigner 
function of a mixed state. This effect, should also be present in the 4-D phase space, but it is orthogonal to the ring 
pattern and was not studdied in j|. 

The characteristic hyperbolic pattern of the mixed Wigner function in the neighbourhood of unstable periodic points 
is the main object of this paper. We shall restrict the treatment to cat maps, i.e. the quantization of linear symplectic 
mappings of the torus. Even though these have nongeneric spectral properties for a chaotic system pX| ] , the eigenstates 
themselves are generically ergodic || in the sense of Von Neumann. Moreover, it was shown by Keating |l2[ that in 
' this case the semiclassical theory is exact, making these maps an ideal probe into the way that the hyperbolic pattern 
emerges as the eigenvalue window that defines the mixed state is augmented. In the limit where we focus on a single 
state, the interference from many patterns washes out any visible structure, but for a sufficiently mixed state, the 
individual pattern is clearly discernible even when the window is not centered on a Bohr-quantized eigenvalue. 

Saraceno and Voros have already shown how the stable and unstable manifolds of the periodic points of the baker's 
map stamp their mark on the various iterates of the corresponding quantum propagator |Q. For this end, they 
employed the Kirkwood phase space representation, which is more appropriate than the Wigner function in the case 
of the baker's map. As for scars in phase space of individual states or states in a narrow energy window, most studies 
have adopted the Husimi representation [|| , which may be interpreted as a Gaussian smoothing of the Wigner function. 
Since the effective window that convolutes the Wigner function is usually chosen to have linear dimension h 1 ^ 2 , it 
dampens the fine structure extending along the stable and unstable manifolds that we here seek to highlight. Indeed 
we show that the spectral Husimi function reproduces the fringes of the spectral Wigner function, but the amplitude 
is modulated by a Gaussian centred on the periodic point. Though this washes out the hyperbolic structure, the 
absence of long range interference allows more periodic orbits to be detected. 
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In section II we review Keating's definition of the spectral Wigner function for maps [ p^| , introducing the Cayley 
parameterization of the linearized classical motion as in Jll| . We then derive the local hyperbolic form of the semi- 
classical contribution of the neighbourhood of a single periodic orbit to the Wigner and Husimi spectral functions. In 
section III we study the particular case of the cat map, following references [[ll 13 . Not only is the basic semiclassical 



theory exact in this case, but our local linearized approximation is then valid throughout the torus. We can then 
clearly visualize the spectral Wigner function as a superposition of hyperbolic fringe patters. As the phase window 
defining the spectral Wigner function is decreased, more fringe patterns are superposed, but it is only as we come close 
to isolating the Wigner function for an individual eigenstate that the interference of the different patterns impedes 
their visual recognition. This interference is dampened in the Husimi spectral function so that more local scars are 
visible. 

To conclude this introduction, we emphasize that the fringe patterns which wc here analyze lie entirely within the 
scope of standard semiclassical theory, restricted to fairly short times. The interesting questions concerning homoclinic 
reccurences refer to individual states depending on the dynamics for long times Likewise, there is no need to 
attempt orbit ressummations [[7 L These have also been presented for the cat maps even though the standard 
semiclassical theory of Keating [jjj is already finite in this particular case. 



II. SPECTRAL WIGNER FUNCTION 

A quantum map is characterized by its time step unitary propagator U. In terms of the eigenangles (or quasi- 
energies) 9 n of the propagator for the eigenstates \ip n >, the unitary propagator for I steps is 

U l = J2^ en \^n><^n\. (1) 
n 

The Weyl representation of U is a function in phase space, 

U l (x)=J2e iWn Mx), (2) 

71 

where, x = (p, q) and w n (x) is the Wigner function for the eigenstate \ip n >. We now construct the spectral Wigner 
function as 

oo 

W(x,9,e) = U l (x)e ae e- We , (3) 



L=— oo 



so that, using rt2h 



W{x,9,e)=J2w n (x) e il ^- e) e-W £ . (4) 



Note that 
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1 — 2e _e cos 

l= — oo k— — oo 



is a periodic comb of broadened delta functions. A slightly different function was used by Keating 12 by taking in 
(||) a sum in I running from to oo. Our definition naturally generalizes the Poisson sum formula that is recovered 
in the limit e — ► 0. Alternatively, for e » 1, J e (j|) « 1. 
Combining (0) and (||), we then obtain 

W(x,9,e) =J2w n (x)S € (9-6 n ), (6) 

n 

where the sum in k has been removed by taking only angles in the range [0, 2ir]. In the limit e — > the spectral Wigner 
function is a comb of delta functions on the eigenangles n , whose amplitudes are the corresponding individual Wigner 
functions w n {x). For values of e larger than the mean level spacing, several eigenstates contribute to the spectral 
Wigner function in a Lorcntzian-likc smoothing of energy width e. Evidently, the spectral Wigner function is only a 
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particular representation of the density operator combining eigenstates within a narrow range of eigenangles. Thus, 
if h n (x) are the Husimi functions for the individual eigenstates f| , the spectral Husimi function corresponding to @ 
will be 



H(x,6,e) = Y,K{x)8 t {e-6 n ). 



(7) 



Note that, contrary to the spectral Wigner function, H{x 1 9, e) is positive definite, because this property holds for 
each individual Husimi function and S e (9) is a positive function. 

To obtain the semiclassical approximation for the spectral Wigner function, we introduce the semiclassical expression 
for the Weyl propagator, 



det(Mj + 1) 



1/2 



exp 



(8) 



where the sum is taken over all the classical orbits j whose center lies on the point x [||. Then Mj stands for the 

monodromy matrix of the Z'th iteration of the map and S,(x) is the center action for the orbit. In the equivalent 
formula for a continuous time system, the phase a — nit, with n = for short times, in which case the sum reduces 
to a single orbit. Generally, the phases for discrete maps must be derived specifically. 
Inserting (||) into (||), we obtain 



W sc (x,9 : e)=Y / £ 



j l=-oo 



|det(M' + 1) 



1/2 



exp 



ue - 
e e 



(9) 



So as to study this object in the neighbourhood of a periodic point xo of the map, we define X = x — xq . The center 
action for the orbit that is close to the periodic orbit is then given by 



Sj (x) = ISq + XBiX + <3(X 3 ), 



(10) 



where ISq = S\(x<q) is the action of the periodic orbit (fixed point). For the choice of 9 = Sq/Ti mod(2ir), we obtain 
the phase coherence of @ at the periodic orbit, which is said to be Bohr-quantized at this angle. Bi is the symmetric 
matrix such that 



with 



Thus Bi is the Cayley parameterization of Mj. 
We can see that the term 



JB X = 



J = 
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-Mj 
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+ Mj 
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|det(Mj 



1 1/2 



2 cosh 



(11) 



(12) 



(13) 



where A is the stability exponent of the orbit. Thus the eigenvalues of the symplectic matrix Mj are exp(— A) and 
exp(A), corresponding to the stable and unstable directions respectively. Using these directions as coordinate axes, 
we obtain 



JB l = 



tanh 








- tanh (^) 



(14) 



so that the contribution to the semiclassical spectral Wigner function from the orbit that lies close to the periodic 
orbit takes the form 



W J sc (X,9,e) = 2J?e e 1. 



- exp [U(^-9) - fpVtanh^f)] ^ _ 1 
f-' 2 cosh (M) ' 4 



(15) 
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where X — (p',q'), the coordinates along the stable and unstable directions respectively. This last expression shows 
that the phase coherence along the stable and unstable directions also holds along each successive hyperbola that has 
the stable and unstable manifolds as asymptotes. This is only a local approximation for a general nonlinear map, 
but it is exact for the cat map in section III. To perform the summation in (|l5|), we approximate tanh = 1 an d 



2 cosh ( 



exp (4m for large values of the XI and I > 0. This leads to a geometrical series, so that 



W 1 sc (X,0,e) = 2Reie l 



cxp(-fpY) 



oo 



1 — exp \i 



'So 

. h 



-0 + i(e+|))] 



2=0 



where the corrections, 



C l (X,0,e)=exp 



il(^r-9)-h 



' exp (- Zip'q' tanh ( ^) ) 
2 cosh (4/) 



/ XI 2i . , 



(16) 



(17) 



decrease rapidly as I grows. 

The dependence of the spectral Wigner function on the phase space variables is reduced to the product p'q' , i.e. it 
is constant along the hyperbolae that contain the classical trajectories. The dominant first term in the bracket of ( |l6| ) 
specifies phase oscillations away from the asymptotes, i.e. the stable and unstable manifolds. The amplitude of these 
oscillations decreases with increasing A and will be maximal for corresponding to the Bohr-quantized orbit. We shall 
show in section IV that the hyperbolic structure predicted by ( |l6| ) is clearly discernible in spite of the interference 
due to other orbits that contribute to the spectral Wigner function. 

In contrast, the Husimi spectral function, 



H(x,6,e) 



dx'W(x',9,e)w x (x'), 



smooths over the spectral Wigner function W(x' , 8, e) with the Gaussian window 

w x {x') = J-exp {~(q> qf ±tf -p) 2 ), 



(18) 



(19) 



which is itself the Wigner function for the coherent state that has (p, q) as its average momentum and position || . 
Generally, the axes of the ellipses in the Gaussian need not coincide with the stable and unstable directions for a 
periodic point. However, in this simplest case we can immediatly evaluate the integral (p"8[). Noticing that each term 
of (|l^) is proportional to 



Va(X) 



- 2 7TPV 



(20) 



with a is 1 or tanh (^)- We obtain each corresponding term of the spectral Husimi function as proportional to 



H a (X) = 



1 + a 2 



exp 



h 1 



jh l 



TP 



Re < exp 



h l 



(21) 



where a = 1 for the dominant resummed term. Thus, in this case, the Husimi spectral function basically dampens the 
hyperbolic Wigner pattern with a Gaussian amplitude modulation centered on each periodic point. The presence of 
the stable and unstable directions is suggested, but the scar is localized in the neighbourhood of each periodic point. 
Qualitatively, we expect the same behaviour in the general case where the axes of the hyperbolae do not coincide 
with the (p',q') coordinates that define the coherent state basis. Though the interference fringes still exist, they will 
be severely dampened far from the periodic orbit. This may allow for the detection of a greater number of periodic 
orbit scars. 



III. THE CAT MAP 



We now apply the present theory to the cat map i.e. the linear automorphism on the 2-torus generated by the 2x2 
symplectic matrix A4, that takes a point a;_ to a point x+ : x + = M,X- mod(l). In other words, there exists an 
integer 2-dimensional vector m such that x + = M.X- — m. Equivalently, the map can also be studied in terms of the 
center generating function |l3| . This is defined in terms of center points 
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X+ + X- 



(22) 



and chords 



where 



£ = x + ~ x - 



-J 



dS(x, m) 
dx 



S(x, m) = xBx + x(B - J)m + ^m(B + J)m 



(23) 



(24) 



is the center generating function. Here B is a symmetric matrix (the Cayley parameterization of M., as in (|lT|)), while 

J = ■ (25) 



We will study here the cat map with the symplectic matrix 



M = 



2 3 
1 2 



, and symmetric matrix B = 



-1 
k 



(26) 



This map is known to be chaotic, (ergodic and mixing) as all its periodic orbits are hyperbolic. The periodic points 
Xi of integer period I are labeled by the winding numbers m, so that 



xi 



1 



(27) 



The first periodic points of the map are the fixed points at (0,0) and (i, i) and the periodic orbits of period 2 are 
[(0,i) , (i,0)], [(±,±) , (§,§)], [(0,±) , (i,|)], [(0,|), (±,f)] and [(0, |), (0, §)]. The eigenvalues of M are e~ A and 
e A with A = ln(2 + \/3) w 1.317. This is then the stability exponent for the fixed points, whereas the exponents must 
be doubled for orbits of period 2. All the eigenvectors have directions v\ = {— \/3, 1) and v 2 — (\/3, 1) corresponding 
to the stable and unstable directions respectively. If we neglect the corrections in (|l^), the semiclassical spectral 
Wigner function has identical fringe patterns for all periodic orbits, except for their overall phase and the amplitude 
of the oscillations, where we must substitute A by nX for an orbit of period n. 

Quantum mechanics on the torus, implies a finite Hilbert space of dimension = — j^-, and that positions and 
momenta are defined to have discrete values in a lattice of separation . The cat map was originally quantized 

by Hannay and Berry [|10[ in the coordinate representation: 



TV 6XP 



— (k 2 



Jk + J 2 ) 



(28) 



where the states < q\qj > are periodic combs of Dirac delta distributions at positions q = j/Nmod(l), with j integer 
in [0, N — 1]. In the Weyl representation |pT| , the quantum map has been obtained in jl3) as 



\det(M + 1)| 



\ ^ e i2TrN[S(x,m)] 



l / ^ 



V2mN \xBx+x(B- J)m+lm(B+J)m] 



|det(7W + 1)| ! 



(29) 



where the center points are represented by x = (-^, j?) where a and b are integers in [0, N — 1] for odd values of N 
jll]]. There exists an alternative definition of the torus Wigner function which also holds for even N. However, it 
is constructed on the quarter torus and this compactification scrambles the hyperbolic patterns as discussed in the 
Appendix. 

The fact that the M. matrix has equal diagonal elements implies that the B matrix has no off-diagonal elements. 
This property will be valid for all the powers of the map and, using (|2^), we can see that it implies in the quantal 
symmetry 



V l M (p,q) = {V l M (-p,q))* = (U l M (p,-q)Y 



(30) 
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It has been shown |To| that the unitary propagator is periodic (nilpotent) in the sense that, for any value of TV 

■fc(A 
M 

1(2to7t + (f>~ 



there is an integer k(N) such that U^jj = e 1 ^ . Hence the eigenvalues of the map lie on the k(N) possible sites 



exp 



k(N) 



, l<m<k(N). 



(31) 



For k(N) < N there are degeneracies and the spectrum does not behave as expected for chaotic quantum systems. 
In spite of the peculiarities in this spectra of quantum cat maps, it is likely that non-degenerate states are typical of 
chaotic maps, such as very weak nonlinear perturbations of cat maps that are known to have nondegenerate spectra 
p5[ . Eckhardt |L6| has argued that typically the eigenfunctions of cat maps are random. 
Keating [12] has shown that the periodicity of the unitary propagator leads to 



JY 



w m (x)Se 



k(N) 

^)E tf (^*, here 



otherwise 



(32) 



m=l v ' (=0 

The symmetry property ([?0]) conjugated with (|32| ) is transmitted to the Wigner functions for the eigenstates: 

ifn(p, q) = w n (p, -q) = w n (-p, q). (33) 

The Husimi representation on the torus depends on the definition of the periodic coherent state |l7|,[l^], with < p >= P 
and < q >= Q: 



oo , 

< P,Q\q k >= ex P) 2lrN 

j=-oo 



k/N) 2 



2uj 2 



iP{j + Q- k/N) 



(34) 



Then the Husimi function for the state \ip n > is defined as 

h n {x) = Tr{\i> n >< Vv>|P, Q >< P, Q\) 



(35) 



where 



Af-l 



W(P,Q){X) = < <l2b-k\P,Q>< P,Q\<lk > 



-i^2(b-k)a 



(36) 



fe=0 



is the Wigner function for the coherent state. This has been studdied numerically by Bianucci et al. |T^] in the case 
of even TV. In the appendix we show that 



oo oo 



j — — oo k— — oo 



x exp < 2ttN 



N 2 



^[P-i^-^) -^.[Q-T-r-i 



N 2 



(37) 



which can be separated into a term that is a periodic Gaussian in phase space centered on the point X = (P. Q) 
whereas the other terms represent Gaussians centered on x — (X + ™) . By definiting the Wigner function on the 
integer lattice (a and b are integer) none of the four Gaussians on the torus are oscillatory, contrary to the case in [jl8| . 
Thus the average over the Wigner function combines the region near X with those near (X + ™). Fortunately the cat 
map that we chose has periodic points of identical structure with this same displacement. Otherwise the hyperbolic 
structure would be erased even more in the Husimi function. 

Figure 1 shows the spectral Wigner function for a cat map with N = 223. In figure 1(a) we see the single Wigner 
function for the pure state with k = 113. For the figures 1 (b) (c) and (d) we show the spectral Wigner function 
W(x, 9, e) for the angle 9 corresponding to the state with k = 113 but for different quasi-energy windows with 
respectively e = 4A,20A and 40A, where A = 2tt/N is the mean level spacing of the system. These figures show 
clearly the imprints of the stable and unstable manifolds that are plotted for clarity in figure 1(d). For the two largest 



G 



values of e we can clearly see scars of the fixed points and also of the periodic orbit [(0, |) , (5,0)] while for the case 
of smaller e = 4A we can also appreciate scars of the periodic orbit, [(|, i) , (^, |)]. 

In figure 2 we compare the exact spectral Wigner function for a cat map with N — 223 in (a) with the semiclassical 
approximation W 3 sc {x, 8, e), correspondingly for h = l/(2nN), in (d), near the fixed point at (1/2, 1/2) whose action 
is So = |. Figure 2.(b) and 2.(c) show respectively an horizontal and and vertical sections of the objects plotted 
in figures 2. (a) and 2.(b). As we can see, the agreement between the exact dynamical system and the semiclassical 
approximation is very good. The modulation arrising in the exact spectral Wigner function is due to the interference 
with other periodic orbits. Figure 3 corresponds to the Wigner and Husimi spectral functions with the same energy 
windows. By erasing the long range extension of the hyperbolic fringes, it becomes possible to detect period-2 points 
of the period-2 orbit, which are not discernible in the Wigner function. 



IV. DISCUSSION 



It has become established lore that the Husimi representation is the ideal instrument to uncover the underlying 
classical phase space structure coded into a bunch of closely associated quantum eigenstates. However, the Gaussian 
smearing of the Wigner function that generates the Husimi function washes out the characteristic hyperbolic patterns 
that stand out in our computations for the cat maps. Indeed the half- width of the gaussian in our figures corresponds 
to y/~N sa 15 points of the Wigner lattice. The imprint of the classical hyperbolicity of the Wigner function is so clear 
that we can detect it even for quasi-energies that do not correspond to a Bohr-quantized periodic orbit. 

The general features exhibited by our calculations should also be discernible for nonlinear systems as, indeed, our 
deduction was not restricted to cat maps. Though the theory in section II is only local, we conjecture that distorted 
hyperbolae asymptotic to curved stable and unstable manifolds will bear fringes reaching out from the periodic point. 
In the case of a chaotic Hamiltonian for two degrees of freedom, this pattern should emerge in two-dimensional sections 
cutting the periodic orbit at a point. This plane should be transverse to that of the orbit, where the Wigner function 
exhibits the pattern of concentric rings discussed in Q . Thus the full picture of the spectral Wigner function emerges 
as a sequence of rings in the plane of the orbit composed with the hyperbolic fringes receeding from the stable and 
unstable manifolds. 

It is remarkable that in the case of the cat map, the hyperbolic fringes surrounding the fixed point emerge, even 
for quasi-energy smoothing e of a few level spacings. Actually, the exponencial cut-off with time in our definition of 
the spectral Wigner function affords equal treatment to all periodic orbits in the denominator of (f[6|). The reason 
why only one hyperbolic fringe system is visible is that for an orbit of period n the lyapunov exponent is multiplied 
by n i.e. A — > nX. As e increases, the interference from the longer orbits is further diminished, revealing identical 
hyperbolic patterns along period-2 orbits, as predicted by our semiclassical formula. Of course, increasing e further 
eventually washes out all the structure. In this respect, the Husimi dampening of Wigner oscillations allows for the 
visibility of a greater number of local scars. 
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APPENDIX: WIGNER FUNCTION FOR A TORUS COHERENT STATE 



To evaluate a discrete torus Wigner function at x = (a/N, b/N) where a and b are half integers, we introduce the 
coherent state wave function (B4) in the Wigner transform (B6I) to obtain 



N-l oo 

W (P,Q)( X ) = £ £ exp^TtAT 

n=0 j,j ,= — oo 



■ r., ■ -a .2(6 — n) /_, a 
iP(j - ]') + t— [P 



1 

The substitution j' = j" — j now leads to 

N-l oo 



f-Q 



n \ 2 



j + Q 



2b 

N 



W(p,q){ x )=Y1 £ exp^iV 

n=0 j,j ff =—oo \ 
2 



-Mi"-""" v " 



N 



^] - 2 Q(n-Nj) + j" 
N I N J> J 



n + Nj 
N 







2 










('-£) 






~ TV 





F )(n + A r i)+/Pj }, 
so that we can now combine the sum over n and j into a single sum over m = n + iVj 

oo oc ^ 



-oo m— — oo 



where we also take j" — > j and 



JL (L-i 

lo 2 \N 2 



and 



— 9 



Q 2 



-0^ + 2jQ + . 7 2 

















nJ_ 






' b 




— i 








2 N 





We now take the Poisson transform of the sum over m so that 



oo oo „oo r 

( X ) = £ / dx exp i 2nN 



IF, 



j=—oo h— — oo 
oo oo 



-i 2 + f Ba + — ) a; - Cj 



j=—oo k— — oo 



x exp < 27tA^ 



N 



N 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 
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Thus we obtain the superposition of four periodic gaussians, one centered at (P, Q) and the other three displaced 
by half integer vectors as described by Bianucci et al [fj"8f . However, it is now necessary to distinguish alternative 
realizations of the Wigner function. In the case of even N, as used in jL8|, the Wigner function is defined on the 
quarter torus (ll), so this is the range of the half integer vectors x\ in the definition of the Husimi function (^5|). By 
centering the quarter torus near the point (P, Q) the Husimi function ( |35| ) is then dominated by the non-oscillatory 
Gaussian, whereas it is only the tails of the other three Gaussians that penetrates in the chosen quarter. 

We have here followed the alternative course of working with odd N which enables us to define the Wigner function 
on an integer lattice on the whole torus. This leads to intelligible scar patterns that would be folded into the quarter 
torus by the rule |ll|] 

W (x + = {-l)^ a+2kh+ ' jkN W{x), (44) 

which would spoil phase coherence along the hyperbolae. 

In the case of the Wigner function for the coherent state, the phase factor in ([5]), which coincides with that in 
(p3|), leads to four non-oscillatory Gaussians, supported by the integer lattice. The Gaussians centered on (P,Q), 
(P + 5, Q) and (P, Q + \) are positive, whereas (P + s, Q + |) centres a negative Gaussian. All four Gaussians must 
then be taken into account in evaluating (|3j]) and it is only the periodic repetitions of these Gaussians whose tails 
can be neglected. 
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FIG. 1. Wigner functions for the state k = 113 with N = 223. 

(a) individual eigenfunction, 

(b) Spectral Wigner Function W(x, 9,e) for the angle 9 corresponding to the quasi-energy of the state with k = 113 and e = 4A, where 
A = 2tt/N is the mean level spacing of the system. 

(c) Idem (b) but with e = 20A, 

(d) Idem (b) but with e = 40A, the black strait lines indicate the stable and unstable manyfolds. 
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(a) (b) 




FIG. 2. detail of the spectral Wigner function with N = 223 near the fixed point [(i, i)] for an energy window of width e = 40A 
centered at the energy of the state k = 113. 

(a) Exact result for the cat map. 

(d) Semiclassical approximation Wg C (x,0,e). 

(b) Wc compare sections of the exact and the semiclassical Wigner functions for q = 0.5 (horizontal section). 

(c) Idem (b) but for p = 0.5 (vertical section). 
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FIG. 3. (a) Spectral Wigncr functions for an energy window centered at the state k = 113 which is not Bohr-quantized with e = 20A 
(b) Corresponding Spectral Husimi Function. The relative maxima along the central vertical line correpond to periodic orbits of period 
two. 
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